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ABSTRACT 

The limits on the mass of the black hole in 23 BL Lac objects are obtained from their luminosities of the broad 
emission line H/3 on the assumption that broad emission lines are emitted from clouds ionized by the radiation of 
the accretion disk surrounding a black hole. The distribution of line luminosity L-ap of all these BL Lac objects 
suggests a bimodal nature, although this cannot be statistically proven on the basis of the present, rather small 
sample. We found that standard thin disks are probably in the sources with Lh/3 > 10"*' erg s"'. The central black 
holes in these sources have masses of IQ^'^^Mq, if the matter is accreting at the rate of 0.025MEdd- For the sources 
with Lh/3 < lO'*^ erg s"', the accretion flows have transited from standard thin disk type to the ADAF type. The 
lower limits on the mass of the black hole in these sources are in the range of L66-24.5 x 10^ Mq. The results 
support the evolutionary sequence of blazars: FSRQ^LBL^HBL. 

Subject headings: BL Lacertae objects: general — galaxies: active — accretion, accretion disks 
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\. INTRODUCTION 

The masses of the central black holes in quasars are crucial 
in understanding the evolution of quasars. Some different ap- 
proaches are proposed to estimate the masses of the black holes 
in active galactic nuclei (AGNs), such as the gas kinematics 
near a hole (see Ho & Kormendy (2000) for a review and ref- 
erences therein). The black hole mass derived from the direct 
measurements on the gases moving near the hole is reliable, 
but it is only available for very few AGNs. For most AGNs, 
the central black hole masses can be inferred from the veloci- 
ties of the clouds in broad line regions (BLRs) and the sizes of 
BLRs (Dibai 1981; Wandel, Peterson, & Malkan 1999; Kaspi 
et al. 2000). Most BL Lac objects have featureless optical and 
ultraviolet continuum spectra, and only a small fraction of BL 
Lac objects show very weak broad emission lines (Veron-Cetty 
& Veron 2000). It is therefore difficult to estimate their cen- 
tral black hole mass from the kinematics of their BLRs. The 
discovery of a tight correlation between stellar velocity disper- 
sion and black hole mass in nearby galaxy bulges offers a new 
tool for estimate of the black hole mass (Ferrarese & Merritt 
2000; Gebhardt et al. 2000). The black hole mass can be esti- 
mated from the measurement of bulge velocity dispersion us- 
ing the M-a relation. The masses of the black holes in some 
nearby BL Lac objects have been estimated from the measure- 
ment of their bulge velocity dispersion (Falomo, Kotilainen, & 
Treves 2001, 2002; Barth, Ho, & Sargent 2002). Wu, Liu, & 
Zhang (2002) adopted the fundamental plane for ellipticals to 
estimate the velocity dispersion and then the hole masses for 
some AGNs. Fan, Xie, & Bacon (1999) derived the masses of 
the black holes in Mkn 501 and some other blazars from the 
gamma-ray variability time scale. They found that the masses 
are around 1O^M0. 

Ghisellini et al. (1998) used a large sample of blazar broad- 
band spectra to study the blazar sequence. They suggested a se- 
quence: HBL^LBL^FSRQ. This sequence represents an in- 
creasing energy density of the external radiation field that leads 
to an increasing amount of Compton cooling. The decrease 
of the maximum energy in the electron distribution causes 
the synchrotron and Compton peaks to shift to lower frequen- 



cies. Georganopoulos, Kirk, & Mastichiadis (2001) argued 
that the radiating jet plasma is outside the broad line scatter- 
ing region in weak sources and within it in powerful sources, 
and the model fits to the spectra of several blazars proposed 
a sequence: FSRQ— ^LBL-^HBL. The evolutionary sequence: 
FSRQ^LBL^HBL, has recently been suggested by D'Elia & 
Cavahere (2000) and Cavaliere & D'EHa (2001). In this evolu- 
tionary sequence, less gas is left to fuel the central engine for 
BL Lac objects, and the advection dominated accretion flows 
(ADAFs) may be in most BL Lac objects. The blazar spectral 
calculations support this scenario (Bottcher & Dermer 2002). 

There are several tens of BL Lac objects in which one (or 
more) broad emission line has been detected. It is therefore 
possible to infer the central ionizing luminosity through their 
broad line emission for these BL Lac objects. The limits on the 
central black hole mass can be obtained, if the accretion type 
in the central engine is known. The cosmological parameters 
Hq = 75 kms"' Mpc"' and =0.5 have been adopted in this 
work. 

2. ESTIMATE ON THE IONIZING LUMINOSITY 

It is not possible to measure the ionizing luminosity directly 
from observations on BL Lac objects, since the observed con- 
tinuum emission from the jets is strongly beamed to us. In this 
case, the optical emission line luminosity can be used to es- 
timate the central ionizing luminosity (Rawlings & Saunders 
1991; Celotti, Padovani, & GhiseUini 1997; Cao & Jiang 1999). 

We can estimate the ionizing continuum luminosity La.Ioh at 
the given wavelength Aq as 



iA.ion(Ao) = 



^line 



(1) 



where Ao is the wavelength of the line, EWtoa is the equivalent 
width of the broad emission line corresponding to the ioniz- 
ing continuum emission (different from the observed contin- 
uum emission). 

The uncertainty in Eq. (1) is the value of EW^oa, which may 
be different for individual sources. We estimate the value of 
EWion from the Boroson & Green (1992) sample. This sample 
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contains all 87 PG quasars (z < 0.5) with high quahty optical 
spectra. The average value of EWh/s is 100 A for 70 radio-quiet 
quasars. Unlike blazars, the optical continuum of radio-quiet 
quasar has not been contaminated by the beamed synchrotron 
emission from the jet. We will take EWion = 100 A and use the 
broad emission Une H/3 to estimate the ionizing continuum lu- 
minosity of BL Lac objects. The BL Lac objects seem to follow 
the statistical behavior of quasars in the correlations between 
radio and broad line emission (Celotti, Padovani, & Ghisellini 
1997; Cao & Jiang 1999). It may imply that the properties of 
BLRs in BL Lac objects are not significantly different from that 
in quasars. If the EWioo of BL Lac objects deviates from that 
of quasars systematically, then the estimated black hole mass 
could be modified with EWion (see further discussion in Sect. 
5). 

3. ESTIMATE OF THE BLACK HOLE MASS 

For a standard thin disk, the ionizing continuum luminosity 
is mainly determined by the central black hole mass and accre- 
tion rate. For a low accretion rate, the accretion flow will transit 
from standard thin disk type to the ADAF type (Narayan & Yi 
1995; Yi 1996). 

3.1. Standard thin accretion disks 

The standard thin accretion disks are thought to be in most 
quasars (Koratkar & Blaes 1999). In this case, the luminosity at 
optical wavelength can be related with the disk luminosity by 
Ld ~ 9ALA,ion(5100 A) (Kaspi et al. 2000). The central black 
hole mass Mbh is then estimated by Mbh ^ 10"^**(Ld/;7i)MQ, 
where ni = M/MEdd- The ionizing luminosity LA,ion can be in- 
ferred from Lh/3, and we can finally estimate the black hole 
mass from the broad line luminosity Lh/3, if ot is known. For 
OT = 1, the lower Umit of the black hole mass is available. 

3.2. ADAFs 

The transition of the accretion flow from the thin disk type 
to the ADAF type occurs while m decreases to a value below 
fficrit (Narayan & Yi 1995; Yi 1996). One possible mechanism 
for the transition might be the evaporation of the disk (Meyer 
& Meyer-Hofmeister 1994; Liu et al. 1999; Rozanska & Cz- 
erny 2000; Meyer-Hofmeister & Meyer 2001). Gu & Lu (2000) 
suggested that the transition can be triggered by the thermal in- 
stabiUty of a radiation pressure-dominated standard accretion 
disk. 

The spectrum of an ADAF: L\{Mbit,m,a,(3) can be calcu- 
lated if the parameters Mbh, in, a, and the fraction of the mag- 
netic pressure (3 are specified. The parameter /3 is defined as 
Pm = (1 -/3)Aot- We can calculate the spectra of ADAFs us- 
ing the approach proposed by Mahadevan (1997). For the fixed 
black hole mass Mbh, the optical luminosity L\ increases with 
a. For AGNs, the viscosity a could be as high as 1, as sug- 
gested by Narayan (1996). In this work, We set a = 0.3 to 
calculate optical luminosity L\ (Yi 1996; Choi, Yang, & Yi 
2001). The fraction of magnetic pressure /3 = 0.5 for equipar- 
tition cases. In fact, our numerical results show that the maxi- 
mal optical luminosity L"^ always requires /? = 0.5, if all other 
parameters are fixed. The accretion rate m < m„h — 0.28a^ 
should be satisfied for an ADAF (Mahadevan 1997). So, vary- 
ing the accretion rate m, we can find the maximal optical lu- 
minosity at Ao numerically for a = 0.3 and (3 = 0.5 (see Fig. 
1). We plot the maximal value of XL'^^{Xo) varying with black 



hole mass Mbh in Fig. 2. We can obtain a lower limit on the 
mass of the black hole from Lh/3 using the relation AL"^ - Mbh 
plotted in Fig. 2. 

For BL Lac objects, most accretion power may be carried by 
strong jets and the accretion flows may probably be described 
by the adiabatic inflow-outflow solutions (ADIOSs) (Bland- 
ford, & Begelman 1999). In this case, the flow is fainter than the 
pure ADAF considered here. The maximal AL"^(Ao) is there- 
fore still vaUd for an ADIOS. It will not affect our estimate of 
the black hole mass Umits. 



4. MASSES OF BLACK HOLES IN BL LAC OBJECTS 

Donato et al. (2001) complied a sample of blazars observed 
in the X-ray band, which includes almost all HBLs, IBLs, and 
many LBLs. We collect all BL Lac objects in their sample and 
1 Jy-i-S4-i-S5 catalogues (Stickel, Meisenheimer, & Kiihr 1994; 
Stickel, & Kiihr 1994, 1996) as our start sample. This sample 
includes most identified LBLs and HBLs. We search the lit- 
erature for all sources with broad emission fine fluxes, and this 
leads to a sample of 23 sources (listed in Table 1). Most of them 
are LBLs, except three HBLs: 0651+428, Mkn 421, and Mkn 
501. The sources with relatively strong broad emission lines 
(EW > 10 A) are considered as genuine HPQs by Veron-Cetty 
& Veron (2000) and are therefore not included in our sample. 
The source 3C 279 (1253-055) is classified as a BL Lac ob- 
ject in Veron-Cetty & Veron (2000) due to its small broad line 
equivalent width, though this source is usually regarded as a 
quasar in other literature. 

Many sources in this sample have detected Mg ii or Ha line 
emission instead of H/3. We use line ratios reported by Francis 
et al. (1991) to estimate H/? line luminosity Lh/3 from the line 
luminosity of Ha or Mg II. The broad emission line luminosity 
Lh/3 are listed in Table 1 . The distribution of the line luminosity 
Lnp is plotted in Fig. 3. 

We use Eq. (1) and the line luminosity Lh/3 to estimate the 
ionizing continuum luminosity at 4861 A. The central black 
hole masses Mbh,i and Mbh,2 can be estimated in the cases of 
a thin disk and an ADAF, respectively. The mass derived for 
the standard thin disk case also depends on accretion rate m. 
The lower hmit of Mbh,i can be available by setting m = 1. If 
the transition of the accretion flow from the thin disk type to the 
ADAF type occurs at m ^ nicnt, we can have an upper limit on 
Mbh,i setting m = 0.025 for a = 0.3. For ADAFs, Mbh,2 is the 
lower limit, since the maximal optical continuum luminosity is 
calculated for the given black hole mass. We list the derived 
black hole masses in Columns (6) and (7) of Table 1 . 
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Fig. 1 . — The optical luminosity of ADAFs at 4861 A varies with accretion FiG. 3. — The distribution of the broad emission line luminosity Lh/3. 

rate rh for different black hole masses (a = 0.3 and /3 = 0.5 are adopted). 
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Fig. 2. — The maximal optical luminosity of ADAFs at 4861 A as functions 
of black hole mass for different values a: a =0.1 (dashed line), 0.3 (soUd), and 
1 (dotted). 



5. DISCUSSION 

The distribution of line luminosity Lug of all these BL Lac 
objects suggests a bimodal nature, although this cannot be sta- 
tistically proven on the basis of the present, rather small sample 
(see Fig. 3). We define the sources withLn/s < 10''^ ergs s~^ as 
population A, and all others are in population B. It is still not 
clear whether this bimodal nature of line luminosities would 
hold for a larger sample. It would be interesting to perform 
high-sensitivity optical observations on a large sample of BL 
Lac objects. 

For thin disk cases, the lower limits on the hole mass would 
be in 10"*"** Mq, while the upper limits on the hole mass would 
be in 10*"'° Mq if m ~ 0.025. For accretion rate m < 0.025, 
the accretion flow would be in ADAF state. In this case, the 
lower limits on the black hole mass are in: 10*^-'^ Mq, which 
is similar to the results of Wang, Xue, & Wang (2001). Noting 
that this is the lower limit, the black hole mass could be much 
higher if the accretion rate is low or/and the viscosity a is small. 

The masses of black holes in all sources of population B are 
in 10''"^ Mq if m ^ 1. It would be interesting to pay much at- 
tention on the seven sources in population A. If the accretion 
flows in these sources are in thin disk state, their central black 
holes would have masses 10^* Mq for any value of m. If the 
accretion rate m is lower than the critical value merit, the accre- 
tion flows in these sources are in the ADAF state. The lower 
limits on the mass of the black hole in these sources are in the 
range of 1.66-24.5 x 10^ Mq. Falomo, Kotilainen, & Treves 
(2001) found that the mass of the black hole in Mkn 501 is 
3.2 X 10^ Mq from the measurement of bulge velocity disper- 
sion, while Barth, Ho, & Sargent (2002) derived the mass of 
(0.9-3.4)x 10^ Mq. Considering fliat our estimate (1.66 x 10^ 
Mq) is the lower limit, the central black hole mass can be higher 
than the limit if the accretion rate is sufficiently low or/and the 
viscosity a is small. Our result is consistent with their results of 
Mkn 501 derived from bulge velocity dispersion. It may proba- 
bly that the sources in population A have already been in ADAF 
state. 
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The sources in population B may have standard thin accre- 
tion disks surrounding the black holes, otherwise some black 
holes should be at least as huge as IO^^Mq. If the accretion 
rate m of these sources is as small as ~ 0.025, slightly higher 
than the critical value below which the accretion flow would be 
in ADAF state, the black hole mass would be in 10^"'° Mq. 
As the fact that most FSRQs have black hole mass higher than 
10^ Mq (Laor 2000; McLure & Dunlop 2001; Gu, Cao, & Jiang 
2001), this is compatible with the unified models of radio-loud 
quasars (Urry & Padovani 1995). 

There are three HBLs in our sample in population A, and all 
sources in population B are LBLs. It may imply that the accre- 
tion flows in all HBLs are in ADAF state. The fact that no broad 
emission line has been detected for most BL Lac objects may 
imply that only a small fraction of LBLs have optically thick 
standard thin accretion disks surrounding the black holes with 
very low accretion rate close to Wcrit. We speculate that these 
LBLs will finally exhaust the gas near the hole and the disks 
will transit to ADAFs. Most other LBLs and HBLs without any 
broad emission Une detected may be in population A, and the 
accretion flows have already been in ADAF state. Otherwise 
the black holes in these sources should be very small, if the 
accretion flows are in standard thin disk state. It is most prob- 
ably that the BL Lac objects studied here with one (or more) 
broad emission Une detected are in the intermediate state of the 
evolutionary sequence from FSRQ to BL Lac object. The fact 
that no HBL is in population B may imply that the evolution- 
ary sequence of BL Lac objects should be LBL— >HBL. The 



results present in this Letter support the evolutionary sequence 
FSRQ^LBL^HBL suggested by D'Elia & Cavaliere (2000). 
If this is the case, then the ratio of the BL Lac objects in pop- 
ulation B to the remainder offers a clue to study the detailed 
evolutionary history of blazars, and it can also be a useful test 
on ADAF models. 

If the evolution of blazars is really regulated by the gas near 
the black hole, the reprocessing optical depth of the BLR would 
decrease with the depletion of the gas near the hole, and the line 
EW\on of the blazar would also decrease along the evolutionary 
sequence. The lower limits on the mass of the hole in evolving 
blazars should be modified with varying EWion- In this work, 
we used a single EWioa to derive the masses of holes in blazars 
and found that the holes have similar masses for these two pop- 
ulations of BL Lac objects. In this case, the lower limits on the 
mass of the hole in the blazars evolving at a later stage would 
become systematically higher than that derived in this work. 
The derived hole masses seem to be consistent with the evolu- 
tionary sequence FSRQ^LBL— >HBL. 
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Table 1 



Broad emission line data and black hole masses 



Source 


Redshift 


Line 


log Liii3 


References 


log Mtb/ 


log Mi,h,2'' 


0235+164 


0.940 


Mgn 


42.2 


CJ99 


8.30 


11.30 


0537-441 


0.896 


Mg II 


43.4 


CJ99 


9.58 


12.60 


0651+428' 


0.126 


Hct 


40.5 


M96 


6.59 


9.39 


0814+425 


0.258 


Mg 11 


40.3 


CJ99 


6.38 


9.15 


08204-225 


0.951 


Mgll 


41.4 


CJ99 


7.51 


10.43 


0823+033 


0.506 


Mgn 


41.8 


CJ99 


7.90 


10.86 


0851+202 


0.306 


H/3 


41.3 


CJ99 


7.44 


10.35 


0954+658 


0.367 


Ha 


41.0 


CJ99 


7.16 


10.03 


1101+384' 


0.031 


Ha 


39.9 


COO 


6.05 


8.78 


1144-379 


1.048 


Mgll 


42.7 


CJ99 


8.80 


11.85 


1253-055 


0.536 


H/3 


42.9 


CJ99 


9.07 


12.14 


1400+162 


0.244 


H/3 


40.5 


SJ85 


6.59 


9.39 


1538+149 


0.605 


Mg 11 


41.5 


CJ99 


7.60 


10.52 


1652+398' 


0.0337 


Ha 


39.4 


COO 


5.57 


8.22 


1722+401 


1.049 


Mgll 


41.5 


WOO 


7.63 


10.55 


1749+096 


0.320 


Mgll 


42.2 


WOO 


8.28 


11.28 


1803+784 


0.684 


H/3 


42.9 


CJ99 


9.99 


12.06 


1807+698 


0.051 


H/3 


39.9 


SJ85 


6.03 


8.74 


1823+568 


0.664 


H/3 


41.7 


CJ99 


7.80 


10.74 


1921-293 


0.352 


H/3 


41.9 


JB91 


8.04 


11.01 


2029+121 


1.215 


Mg II 


42.2 


CJ99 


8.29 


11.29 


2200+420 


0.068 


Ha 


39.5 


SJ85 


5.63 


8.28 


2240-260 


0.774 


Mgll 


41.8 


CJ99 


7.98 


10.94 



'in unit of 0.025Mq /m 
''in unit of Mq 

'hbl 



References. — COO: Cao (2000); CJ99: Cao & Jiang (1999); JB91 : Jackson & Browne (1991); 
M96; Marcha et al. (1996); SJ85: Sitko & Junkkarinen (1985); WOO; Veron-Cetty & Veron 
(2000). 



